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We study the effects of single-impurity scattering on the local density of states in the high- Tc 
cuprates. We compare the quasiparticle interference patterns in three different ordered states: d- 
wave superconductor (DSC), d-density wave (DDW), and coexisting DSC and DDW (DSC-DDW). 
In the coexisting state, at energies below the DSC gap, the patterns are almost identical to those in 
the pure DSC state with the same DSC gap. However, they are significantly different for energies 
greater than or equal to the DSC gap. This transition at an energy around the DSC gap can be used 
to test the nature of the superconducting state of the underdoped cuprates by scanning tunneling 
microscopy. Furthermore, we note that in the DDW state the effect of the coherence factors is 
stronger than in the DSC state. The new features arising due to DDW ordering are discussed. 
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I. INTRODUCTION 

One of the intriguing features of the cuprate high tem- 
perature superconductors is the existence of a pseudo- 
gap phase in the normal statei Elucidating the nature 
of the pseudogap would be an important step in under- 
standing the physics of high-temperature superconduc- 
tivity. Among many scenarios, one concrete proposal is 
that the pseudogap is due to a hidden broken symmetry^ 
of d^2_y2 type in the particle-hole channel.'^ In this paper 
we examine the consequences of the presence of DDW or- 
dering on the quasiparticle interference patterns observed 
by scanning tunneling measurements (STM). 

Recently, interesting STM measurements^*^ on Bi- 
2122 have been performed to obtain Fourier transformed 
spectra of the local density of states. This technique is 
called Fourier transform scanning tunneling spectroscopy 
(FT-STS). So far, the experiments have been carried 
out in the superconducting state at low temperatures, 
and the results can be interpreted as interference pat- 
terns due to elastic scattering of the quasiparticles from 
impurities Additional support of this interpretation 
comes from a more recent work in Ref. The scatter- 
ing of quasiparticles between regions of the Brillouin zone 
with high densities of states yields peaks in FT-STS. In 
a c?-wave superconductor (DSC), the regions of high lo- 
cal density of states (LDOS) are situated at the tips of 
the banana-shaped contours of the quasiparticle excita- 
tion spectrumi4 The wavevectors of the observed peaks 
are consistent with the wavevectors that connect the tips 
of these bananas (see Fig^J. Thus, the experimental re- 
sults are consistent with the known pictures of the Fermi 
surface and the LDOS in the pure DSC state. Here, 
we investigate the possible form of STM spectra in both 
the superconducting and normal states from the perspec- 
tive of DDW. The T-matrix formalism)^ combined with 
a numerical analysis, are used to obtain the quasiparticle 
interference patterns of the various ordered states. 

We first consider a pure DSC state and recover results 
similar to Ref. 0. Though the band structure we used 



is different from theirs, the main features such as the 
emergence of peaks at particular wave vectors are still 
observed. However, as we will discuss later, differences 
in the band structure can have large effects on the details 
of the STM spectra, especially due to coherence factors. 
Fortunately, none of our robust qualitative conclusions 
depend on the spurious sensitivity to the band structure. 

Our main focus is to study the mixed, or coexisting, 
DSC-DDW and the pure DDW states. According to 
the DDW theory,^ the mixed DSC-DDW and the pure 
DDW characterize the underdoped superconducting and 
the pseudogap phase respectively. We first analyze the 
mixed state and find that if we probe the system at en- 
ergies larger than the DSC gap, significant differences 
among the interference patterns emerge compared to a 
pure DSC state, while at energies below the DSC gap, 
the patterns are almost identical. For example, let us fo- 
cus on the position of the peaks along the (tt, tt) direction 
in the interference patterns. For a pure DSC state they 
should move away from the origin with increasing energy 
and eventually get close to the (tt, tt) point for energies of 
the order of the gap. If the state is a mixed state (DSC- 
DDW) with two different order parameters, the total gap 
is roughly equal to the square root of the sum of the two 
gaps square. In this case the same peaks should shift 
away from the origin faster (for lower energy) and get 
close to the (tt, tt) point at energies of the order of the 
DSC gap which is typically much smaller than the total 
gap. For energies larger than the DSC gap, the main fea- 
tures of the STM spectra is significantly different from 
those in the pure DSC state. Therefore, the energy of 
the DSC gap in the mixed states marks a transition in 
the interference patterns. The effect becomes less pro- 
nounced with an increase in doping, as the DSC gap is 
getting closer to the total gap of the system. Based on 
this observation we believe that the FT-STS obtained 
below the superconducting transition temperature in the 
underdoped cuprates may reveal important information 
about the nature of the order parameter. 

In the pure DDW state, the interference patterns re- 
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veal features similar to the ones observed in the DSC 
state and others that are entirely new. These features 
could be used to detect DDW order. One of the impor- 
tant differences from the case of DSC order is that the 
DOS is not as strongly peaked in the DDW state, so the 
STM spectrum is more strongly dependent on coherence 
factors. These depend not only on the order parame- 
ter and on the band structure, but also on the details of 
the scattering and on the type of impurities. We illus- 
trate this effect by comparing the interference patterns 
for two parameter sets which are identical save for the 
sizes of their DDW gaps. Changing the size of the DDW 
gap has relatively little effect on the equal energy con- 
tour plots. However, the interference patterns are dra- 
matically different due to the difference in the coherence 
factors. The authors of Ref. emphasized their im- 
portance. We show that their observation is pertinent 
and, in fact, explains why their lowest-order Born ap- 
proximation results miss some of the physics captured 



in the full T-matrix approach. Also, the strong depen- 
dence of the STM spectra on coherence factors suggests 
that they could be used to extract detailed information 
about the DDW gap, band structure, etc., which could 
be compared with angle-resolved photoemission spectra 
(ARPES) data. 

The paper is organized as follows. In section II we 
describe our model and the analytic formalism. In section 
III, we describe the numerical analysis and discuss the 
results; we conclude in section IV. 



II. ANALYTIC FORMALISM 

In the presence of DDW and DSC order, the general 
mean field Hamiltonian for a high superconductor is 
given by 
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where W{k) — Wo(cos/ca; — coaky)/2 and A(fc) = 
Ao(cos kx — cos ky)/2 are the DDW order (see also ref. j?} 
and DSC order parameters respectively, and Q = (tTjTt). 
The sums over k include all the wave vectors in the first 
Brillouin zone (BZ), jfe^,! < tt, \ky\ < n. We take the 
lattice constant a = 1. 

The kinetic energy is e(fc) = ei(fc) + £2(^)1 where ei/e2 
are respectively odd/even with respect to the shift of the 
wavevector by Q; ei(fc) = —ei{k + Q) and e2{k) = e2{k + 
Q). In the numerical calculation, we focus on a simple 
t,t' band structure with dispersion ei(fc) = t{coskx + 
cos ky)/2 and £2 (fc) = t' cos k^ cos ky . 



J 



To simplify the calculation, we introduce a 
four-component (spinor) field operator, = 

In this new basis, the 



above Hamiltonian can be written as 



H = J2ij+A{k)i;, 



(2) 



where fc is summed over half of the original Brillouin zone 



(reduced Brillouin zone - RBZ), namely, \kx\ 
Ak is a four by four matrix given by 



\ky\ < TT. 



A, 



ei(fc) + e2(fc) iVV{k) A(fc) 

-iW{k) -ei(fc) + e2(fc) -M -A(fc) 

A*(fc) -ei(fc) - e2(fc) +M ^W{k) 

-A(fc)* -iW{k) €i{k) - €2{k) + n 



r 
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The eigenvalues of Ak are ±£'i(fc) and ±i?2(fc), where Up to a constant, we can write it as 
i?i,2(fc) = {(v/ei(fc)2 + W{kf ± [e2(fc) - + A(fc)2}i 

(4) Hin.p= ^kVik,k')^k', 

We consider impurity scattering of the form: 

^imp = Y X! ^kk' afjclaCk' (3 (5) 
k,k' Q,/3=T4 



(6) 



k,k'ert.BZ 



where V{k^ fc') is a four by four matrix. 
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One can define a finite temperature (imaginary time) mulation : 
Green's function, 



G(fci,fc2,T) = -Tr [e-'^(^-")T.Vfc,m+ (0)] (7) 

where K ^ H - fiN, e"'^^^ = Tre"^^, and is the 
imaginary time ordering operator. The impurity scat- 
tering problem can be solved by computing the Fourier 
transform of the Green's function from the T-matrix for- 
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G(fci, fe, icjn) = Goiki,iuJn)T{ki, fc2, iw„)Go(fc2, ia;„),(8) 
where 

Go{k, iLiJ„)~^ = iujrJ - Ak, (9) 
/ is the 4x4 identity matrix, and 



T{ki,k2,iu!n) = V(ki,k2) + ^ V{ki, k')Go{k' ,iu!n)T{k' , k2,iuJn). 



(10) 
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For simplicity, we take the impurity scattering potential 
to be a delta function so that the matrix V is independent 
of k and k'; V{k, k') = V . For this case, we can solve Eq. 
(|1(J|> and obtain 



T{iujn) = [l-V dkGo{k,tij,,)]-W 



(11) 



where the only difference from the standard result is that 
the integral over k is over the RBZ. 

Consequently, the local density of states p{q,uj), is 
given by 



(12) 



kGRBZ 



where g{k,q,uj) is defined as follows. Let k' ^ k + q. If 
k' is in the RBZ, 

4 

g{k, q, uj) = ^[Gii{k, k' , SiUj) - G*i{k' , k, SiUj)] (13) 
1=1 

where Si = 1 for i = 1,2 and = — 1 for i = 3, 4. If k' is 
not in the RBZ, let k" = k + q — Q. For this case 

g{k,q,u;) = ^ [G,,j+i(fc, A:", s^w) - G*j+i(fc", /c, Sjw) 

4=1,3 

+ G,+i4k,k",.s,Lo)-G*+^.,{k",k,s,u)] (14) 

Here G(fci,A:2,w) is obtained by analytical continuation 
iun uj + i6 oi G{ki, k2,iujn) from imaginary frequencies 
to real frequencies. The complexity of the above formula 
stems from translational symmetry breaking in the pres- 
ence of DDW order. 



III. NUMERICAL CALCULATION AND 
DISCUSSION 

We compute the local density of states p{q,uj) using 
the 4x4 impurity scattering matrices. For potential 
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scattering given by a (5-function, they are 

V{k, k') = Vn 
for a non-magnetic impurity, and 
V{k, k') = Vm 



/ 1 1 

110 

0-1-1 

Vo -1 -1 



(15) 




(16) 



for a magnetic impurity. 

The results are reported for the representative values: 
Vn = Vm = 0.1 eV. For the band dispersion, we choose 
t — —1.2 eV and t' = 0.36 eV. The chemical poten- 
tial is selected to be equal to —0.36 eV. The imaginary 
part of the energy S ~ 0.5 meV is used for the entire 
numerical calculation. We have checked that the results 
are unchanged for smaller values of S. Following Ref. |^ 
a 400 X 400 lattice is used in our analysis, and the re- 
sults are displayed in the (— tt, tt) x (— tt, tt) interval on a 
49 X 49 grid for any given frequency. The choice of these 
parameters is representative. We have repeated our cal- 
culations for a number of different set of parameters; the 
conclusions remain unchanged. 



A. The DSC state 

First we analyze the interference patterns in a pure 
superconducting state with Aq ~ 25meV. This has also 
been studied in Ref. P for a different bare band struc- 
ture. The constant energy contour plots are shown in the 
Fig. n^b) where we label the wavevectors, which are ex- 
pected to be associated with the peaks in the interference 
patterns. In Fig. |21 we show a comparison between our 
results and the results in Ref. where a different band 
structure^'^ which is flatter near (tt, 0) is used. Some of 
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the features in our results are the same as in Ref. 0, but 
there are also several clear differences. 

For the case of a nonmagnetic impurity, the peaks as- 
sociated with the wavevectors numbered 3, 4, and 7 are 
observed for a large energy range. However, the intensity 
is high mainly along the diagonal or symmetrically about 
the diagonal, contrary to the results in Ref. Isl where 
regions of high intensity emerge along the (0, ±1) and 
(±1,0) directions when the energy |u;| reaches 15meV. 
The peak associated with the wave vector numbered 1 
appears for some energies, but with weaker intensity than 
the peaks associated with the above three wave vectors. 

For the case of a magnetic impurity, there are fewer fea- 
tures, so it is more difficult to relate the observed inten- 
sity to the scattering from the tips of the banana-shaped 
contours. 



B. The mixed DSC and DDW state 

If DDW order is the origin of the pseudogap phase, 
the mixed state should describe the underdoped cuprates 
at low temperatures. In Fig. ^Jd) we sketch the equal 
energy contours for the band structure of a mixed DDW- 
DSC state. For energies lower than the DSC gap the 
equal energy contours are now four pairs of bananas due 
to the doubling of unit cell by the DDW order, and the 
regions of high DOS are situated at their tips. However, 
for energies larger than the DSC gap, the equal energy 
contours become elliptical, which is characteristic to the 
DDW state. 

In Fig. 121 we plot our results for the mixed state for 
nonmagnetic impurity scattering. We note that for ener- 
gies lower than the DSC gap, we can identify exact fea- 
tures characteristic to the pure DSC state. The patterns 
are almost identical to those obtained in a pure DSC state 
with the same DSC gap Aq. In particular the peaks along 
the (tt, tt) axis disperse with energy in the expected fash- 
ion. However, they reach the corners of the BZ at energy 
comparable to the DSC order par ameter Ap , and not to 
the full gap of the system Ai = \/Wq + Aq. Above the 
energy Aq, a whole range of different features emerge, 
which are characteristic to the DDW state. For compar- 
ison, in Fig. 21 we plot the corresponding results for the 
case of a pure DSC state with a gap equal to the total gap 
of the mixed state Aj. Indeed, the peaks in the pure DSC 
state disperse much slower with energy and can be ob- 
served at higher energies than in the mixed state. Based 
on these observations we expect the FT-STS measure- 
ments to be an experimental tool to observe the coexis- 
tence of DDW order with DSC order in the underdoped 
cuprates. If this state is indeed a coexisting DDW and 
DSC state, measurements done at various energies should 
reveal that the diagonal peaks in the spectra should ap- 
proach the corners of the BZ well before the full gap 
is reached. Also, entirely new features should arise in 
the STM spectra for energies above the energy for which 
the peaks have reached the corners. The new features 



would be similar to what one would expect to see in a 
pure DDW state. Broad regions of high intensity are ex- 
pected rather than sharp peaks, as we will discuss in the 
next subsection. When the doping is increased, since the 
total gap At is roughly constant for different doping lev- 
els, the DDW gap decreases and the DSC gap increases. 
The energy at which the transition described above hap- 
pens in the interference patterns is thus expected to get 
closer and closer to the total gap and eventually equal to 
it when the DDW gap vanishes and the DDW ordering 
disappears. 



C. The DDW state 

We now turn to a pure DDW state. In Fig. ^c) we 
plot a typical band structure of a pure DDW state. We 
note that the equal energy contours are now elliptical and 
we therefore expect the high DOS regions to be situated 
at their tips, though, due to the smaller curvature of 
these contours, the distinction between the high and low 
regions of DOS are not as strong as in the case of DSC. 
In Fig. O we plot the corresponding LDOS results for 
Wo = 25 meV, and Wo = 40 meV. 

Even though one can still track the presence of peaks, 
they are not as pronounced. This may be the result of 
a more uniform DOS. The position of some of the peaks 
(such as some along the (tt, tt) direction and symmetri- 
cally about it, and some along the (7r,0) direction) may 
be traced back to the positions of the tips of the elliptical 
contours in the band structure. These peaks are present 
in most of the pictures, but their intensities vary. 

Also, for the chosen band structure, peaks located 
around (±^,0) and (0, ±^) can be found over a broad 
range of energy. They are present in both the Wo = 
40meV and Wq = 25meV data. In particular, their posi- 
tions hardly change with energy for positive energies. We 
note that their peculiar positions of roughly (±^,0) and 
(0, ±^) is an effect of the particular set of parameters in 
the band structure. The origin of these peaks is the scat- 
tering indicated by the arrow in Fig.^c). Since the equal 
energy contours may change drastically with a change of 
parameters, so may the position of the peaks. We also 
note that, since the equal energy contours change very 
little with energy for positive energies, the peaks also do 
not move when the energy is varied. 

Many other peaks however cannot be explained by sim- 
ple band structure arguments, and their positions change 
drastically from one plot to the next. 

The predominant feature for the DDW state is thus 
the presence of broadly distributed scattering points in 
the interference patterns at low energy. However, lines 
with relatively high intensity are also present in our data. 
Such lines actually occur in all of the states considered. 
Furthermore, the intensity and size of these lines vary 
strongly with even slight changes of parameters. This is 
in contrast to the statement in Ref. [T^ that the presence 
of high intensity lines in the spectrum is a characteristic 
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of the DDW state. However, since the equal energy con- 
tours are elliptical, the scattering patterns in the DDW 
states are much more likely to be affected by coherence 
factors than the scattering patterns in other states. Be- 
cause of this effect, the comparison between DDW and 
DSC in Ref. HI based on Born approximation, should 
not be considered as generic, although we generally agree 
that the interference patterns are different in these states. 
We will discuss this in more detail in the next subsection. 

As a final observation, we note that if we tune the pa- 
rameters such that electron pockets around (±7r, 0) and 
(0, ±7r) appear in the DDW band structure, our calcu- 
lation shows circles of high intensity in the BZ corre- 
sponding to electrons with different momenta scattering 
off each other. This appears to be the dominant fea- 
ture. However, the presence of such electron pockets is 
not consistent with ARPES data-'^'' and we do not include 
the results for this situation in the present work. 

D. Discussion of the coherence factors 

The interference patterns not only depend on the band 
structure and on the order parameters, but also on the 
coherence factors. This dependence on coherence factors 
has quite complex effects and makes the interference pat- 
terns sensitive to many details. 

First, the coherence factors depend on the strength and 
the type of the impurity scattering. In Ref. 9, the in- 
terference patterns from nonmagnetic and magnetic im- 
purities were shown to be so different that entirely dif- 
ferent peaks were observed, even for pure s-wave scat- 
tering. Clearly, similar or more pronounced effects may 
arise in the case of non s-wave scattering. Furthermore, 
the coherence factors depend on the impurity scattering 
strength. This may be implied from the sensitivity of the 
interference patterns to whether we use the Born approx- 
imation or a full T-matrix calculation for a reasonable 
impurity strength. 

Second, the coherence factors are also strongly sensi- 
tive to the band structure and order parameters. The 
band structure dependence causes the dramatic differ- 
ence between our results for the DSC state and the results 
in Ref. 9. The flatness of the band structure of Ref. 0| 
leads to differences in the DOS and also to important dif- 
ferences in the coherence factors. In the DDW state, the 
effect of coherence factors is more dramatic than in the 
DSC states, as seen from the plot of our results for iden- 
tical parameter sets save for the DDW gaps. Though the 
band structures for the two chosen DDW gaps are very 
similar, as seen in Fig.|31there arc considerable differences 
between the corresponding FT-STS patterns. 

Third, the coherence factors are sensitive to energy. 
This has been shown consistently in all our calculations. 
Even for similar band structures, some peaks can only be 
identified for some ranges of energy, while others appear 
or disappear in a manner counterintuitive to expectations 



based on band structure. 

Finally, the coherence factors should also be expected 
to depend on the correlation and distribution of impuri- 
ties in the materials. 

The sensitivity of the coherence factors to so many 
details makes it difficult to extract information about 
the dispersion of quasiparticles from the interference pat- 
terns. It is rather hard to identify a feature that can be 
traced in a broad energy range. This is particularly true 
for the DDW state, where the quasiparticle scattering in- 
terference is much more sensitive to the changes of the 
parameters in the model. 



IV. CONCLUSIONS 

In this paper we have analyzed the effects of scattering 
from a single impurity on FT-STS in a high Tc supercon- 
ductor with DDW ordering. For the case of a pure super- 
conducting state, we recovered results similar to Ref. I^. 
We note the presence of peaks dispersing with energy in 
the expected fashion. However, some important features 
are strongly dependent on the band structure and on 
other parameters in the model even if the banana-shaped 
equal-energy plots are similar. In the case of a mixed 
DDW-DSC state, we note that the DSC order dominates 
the interference patterns at energies lower than the DSC 
gap Aq, while for higher energies, a transition to a spectra 
characteristic of the DDW phase occurs. This is marked 
by a sudden change of the general features of the STM 
spectra around an energy E = Aq. In particular, we 
show that the dispersion of the diagonal peaks in the 
spectrum can be used to identify this transition. In the 
underdoped regime, the DSC gap Aq is reduced when the 
doping decreases, so the energy at which the transition 
happens should also be reduced. Based on this obser- 
vation, we predict that the STM measurements in the 
underdoped cuprates below Tc may reveal important in- 
formation about the nature of the gap. In the pure DDW 
state the results are more complex; the peaks are more 
broadly spaced and not so pronounced and are more sen- 
sitive to changes of parameters. It would be interesting 
to see if STM experiments in the pseudogap regime can 
provide new insight into the problem. 
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FIG. 1: Equal energy contour plots in different states: (a) normal state without any order, (b) d-wave superconducting state 
(DSC), (c) d-density wave state (DDW), and (d) Mixed DSC and DDW state (DSC-DDW). 
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FIG. 2: Fourier transform of STM spectra: (a) in a pure DSC state with the t — t' band structure and Aq = 25 meV, (b) 
in a pure DSC state with the band structure in Ref. ^3 = 25 meV and at energies equal to — 6meV and — 18meV, 

for nonmagnetic impurity scattering {Vn ~ lOOmeV). The differences between (a) and (b) reflect the sensitivity of FT-STS 
spectra to band structure details. 
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FIG. 3: Fourier transform of STM in the coexisting DSC-DDW state with the t — t' band structure, Ao = 25meV and 
Wo = 40meV, for nonmagnetic impurity scattering {Vn = lOOmeV). The peaJis along the directions (±7r, ±7r) are marked by 
circles. As noted, these peaks get close to the corners of the BZ at energies equal to the DSC gap Ao, and not to the total gap, 
which is At w 47meV. Also, the general features of the spectra change dramatically at energies larger than Aq. 
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FIG. 4: Fourier transform STM spectra in a pure DSC state with the t — t' band structure and Ao = 47meV, for nonmagnetic 
impurity scattering Vat = lOOmeV. Again, the peaks along the directions (±7r, ±7r) are marked by circles. In this case, the 
peaks get close to the corners of the BZ for energies equal to the DSC gap Ao = 47meV. 
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FIG. 5: Fourier transform STM spectra in the pure DDW state with the t — t' band structure and Wo = 40mcV, for non- 
magnetic (Vjv = lOOmeV) and magnetic {Vm = lOOmeV) impurity scattering. The peaks at (0, ±7r/4) and (±7r/4, 0) are marked 
by circles and can be observed for a large range of energies. 



